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In order to characterize the behavior of tunable diode lasers (TDL), the first step 
in the project involved the redesign of the TDL system here at the University of 
Tennessee Molecular Systems Laboratory (UTMSL). Having made these changes it 
was next necessary to optimize the new optical system. This involved the fine 
adjustments to the optical components, particularly in the monochromator, to 
minimize the aberrations of coma and astigmatism and to assure that the energy 
from the beam is focused properly on the detector element. The next step involved 
the taking of preliminary data. We were then ready for the analysis of the 
preliminary data. This required the development of computer programs that use 
mathematical techniques to look for signatures of chaos. Commercial programs 
were also employed. We discovered some indication of high dimensional chaos, but 
were hampered by the low sample rate of 200 KSPS (kilosamples/sec) and even 
more by our sample size of 1024 (IK) data points. 

These limitations were expected and we added a high speed data acquisition 
board. We incorporated into the system a computer with a 40 MSPS (million 
samples/sec) data acquisition board. This board can also capture 64K of data points so 
that were then able to perform the more accurate tests for chaos. The results were 
dramatic and compelling — we had demonstrated that the lead salt diode laser had a 
chaotic frequency output. 

' Having identified the chaotic character in our TDL data, we proceeded to stage 
two as outlined in our original proposal. This required the use of an Occasional 
Proportional Feedback (OPF) controller to facilitate the control and stabilization of 
the TDL system output. The controller was designed and fabricated at GSFC and 
debugged in our laboratories. After some trial and error efforts, we achieved chaos 
control of the frequency emissions of the laser. 

The two publications appended to this introduction detail the entire project and 
its results. 

Here we provide the references to the publications: 

Stabilizing Lead-Salt Lasers: Understanding and Controlling Chaotic Frequency 
Emission, Gordon Chin, Larry R. Senesac, William E. Blass, and John J. Hillman, 
Science 274, 1498-1501 (1996) 

Stabilizing Lead-Salt Diode Lasers: Measuring and Controlling Chaotic Frequency 
Emission, Larry R. Senesac, William E. Blass, Gordon Chin, and John J. Hillman, 
Proceedings of the 4th Experimental Chaos Conference, August 1997, Boca Raton, 
FL, in press. World Scientific. 
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Stabilizing Lead-Salt Diode Lasers: 
Understanding and Controlling Chaotic 
Frequency Emission 

Gordon Chin, Larry R. Senesac, William E. Blass, 

John J. Hillman 

Lead-salt tunable diode lasers (TDLs) are the only devices currently available that can 
generate tunable monochromatic radiation at arbitrary wavelengths between 3 and 30 
micrometers and are particularly useful for high-resolution spectroscopy over a wide 
range of spectra! regimes. Detailed observations of TDLs show that the observed in- 
strumental linewidth is actually a temporal average of many narrow (less than 0.5 
megahertz) emission “modes.” The time scale characteristic of these modes, which 
appear to be of relatively constant intensity, is of the order of a microsecond. The laser s 
behavior is highly suggestive of a chaotic process, that is, seemingly random excursions 
of a dynamic variable (frequency) within a bounded range. This report shows experi- 
mentally that TDL emissions are indeed chaotic. Furthermore, in a simple and robust 
fashion, this chaotic behavior has been successfully controlled with the use of recent 
techniques that take advantage of chaos to produce a narrow band laser output. 


Chaotic behavior is an intrinsic property of 
a nonlinear system. Mechanical systems, 
communication and electronic systems, bio- 
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logical systems, the solar system, optical and 
laser systems, and many other nonlinear sys- 
tems often exhibit chaotic behavior. Chaos 
is also often called deterministic random mo- 
tion. The morion is deterministic because its 
trajectory' can be calculated for all times 
given the starting conditions, but it is also 
motion in which instabilities appear every- 
where in the system’s trajectory’ in phase 
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space. rhe coordinate system describing a 
s\Mcm‘> position and velocity. Thoc insta- 
bilities force the trajectories of nearhv phase - 
space points to diverge. This divergence is 
manifested as a sensitivity to initial condi- 
tions, another characteristic of chaotic mo- 
tion. However, the nonlinearities of the sys- 
tem fold the trajectories hack into a confined 
region. In this manner, chaotic motion am- 
plifies small differences, errors, or noise, that 
is, the instabilities stretch phase space, 
whereas the nonlinearities, or folding mo- 
tion, keep the trajectories hounded. Bound- 
ed trajectories signal rhe presence of an at- 
tractor. Another indication of chaotic be- 
havior is the existence of a “strange" attrac- 
tor with its distinctive fractal character. It is 
this fractal characteristic that gives chaotic 
motion its random appearance. 

More recently, investigations of nonlin- 
ear dynamics have moved from purely the- 
oretical and computational studies to exper- 
iment. Some of these experiments show 
that chaotic behavior is a useful feature in 
itself and more importantly may he exploit- 
ed to perform practical functions (/). Fol- 
lowing theoretical studies that reveal meth- 
ods that stabilize unstable orbits (2, 3), a set 
of striking experiments ensued that demon- 
strated that chaos can be controlled in a 
surprisingly simple manner (4, 5) and in as 
diverse a ranee of problems as cardiac ar- 
rhythmia (6), thermal convection loops 
(5), electronic circuitry (7), an optical sys- 
tem (8), and solid-state lasers (9). 

The practical demonstration of chaos 
control led us to consider whether these 
techniques are applicable to tunable diode 
lasers (TDLs). Lead-salt TDLs are valuable 
tools that are used to measure the molecular 
parameters of atmospheric constituents im- 
portant in planetary atmospheres (10). The 
advantages of TDLs are wide spectral cover- 
age (3 to 30 pun), achieved by varying the 
salt composition of an epitaxially grown 
lead-salt crystal (such as Ph,_ x Sn x Te or 
FhSj % Se % ), and high spectral resolution 
( — 1 0 4 cm" 1 ). The TDLs are operated be- 
low 100 K and are tuned by changing the 
temperature or injection current (which ef- 
fectively causes a change in the temperature 
of the laser junction). The combination of 
w ide spectral coverage, continuous tunahil- 
ity. high resolution, high output power, and 
ease of use makes TDLs ideal for the study of 
line positions, line strengths, and pressure 
broadening coefficients, especially of mole- 
cules with closely spaced spectral transitions. 

Despite the extremely successful use of 
TDLs in laboratory spectroscopy, the full 
potential of the TPL has not been realized. 
Careful characterization of a TDL shows 
that the intrinsic hnewidth of its emission is 
v ery narrow (<0.5 MHz or ^ cm 1 ), a 
factor of 10 better than the instrumental 


resolution achieved in labour rv applica- 
tions. Tlie observ ed instrument . unewidrh 
is actually a temporal average : many nar- 
row frequency emissions or “m des” CO, 
1 / ). These instantaneous mode-, which are 
relatively constant in intensity, c.oange on a 
time scale on the order of a microsecond. 
This behavior is a natural consequence of 
the intrinsically nonlinear operation or a 
laser junction, w here many micr 'Copic fac- 
ets or lasing regions may ex or. Minute 
changes in the junction — such o m tem- 
perature, dimension, or index < if retrac- 
tion — can quench one lasing m> -de while 
allowing another to reach threshold. The 
detailed dynamics of such a c triplicated 
system is probably not amenable to analyt- 
ical description. 

Fortunately, numerical me: nods exot 
that can identity the presence r chaotic 
behavior without requiring .* detailed 
knowledge of the underlying dynamics 
However, the challenge rema;:> how to 
capture the trajectory of the TDL tn which 
the dynamical variable, frequency, varies at 
a megahertz rate. This information, is also 
critical in our goal to control me TDL s 
chaotic emission. 

The solution is to construct a rrequenev- 
to-amplitude converter (F/A) wirh a ra>t 
response time. We use one of tw o .ipproach- 



Fig. 1. The frequency variation of a r sc -salt TDL 
operating near 1 029 cm 1 using a lire of ethyiene 
as the F/A convener. Each panel ~a$ 30. COO 
points digitized with a 20-MHz 8-b: snalog-to- 
digital converter. The colored points a r e the mea- 
sured data. The black plot shows re results of 
applying a digit low-pass filter to tre measure- 
ments. The frequency scale was ca mated with 
the use of methods similar to those -$ed in 1 10 . 
7 /}. (A) Chaotic TDL emission withcu: zontrei. B) 
TDL emission with the OPF electro*' :s on 1 C) 
Detector noise level when the \as~' beam is 
blocked. The difference between the * :ered and 
unaltered data is identical to the nose. shovv.rg 
that we succeeded in taking out onlv me detector 
noise. The interval of the measurements was T .5 
ms; the data in tne three panels wee not ta^en 
contiguously. 


l' either .in aKorprion cell v. draining a 
g i 'cou> molecular ^impic or .i rahry- Perot 
ct-.tlon i> placed in the path rhe TDLs 
beam. Bv adjusting the TDL r«> -a>e on the 
'tcep side of either rhe gas line or the 
F ibrv-Peror fringe, changes in rhe TDL fre- 
quency are converted to change- in detect- 
ed >ignal amplitude. The F/A wc used was 
an absorption transition of ethylene (CDH 4 ) 
near 102° cm" 1 (Fig. 1A). However, an 
eta Ion or any other gas line within the 
tuning range of rhe TDL can work just as 
well. The advantages of u>mg a gas absorp- 
tion line are minimal losses in photon flux 
compared with those tor a high-rtnesse eta- 
Ln and rhe ability to control rhe slope of 
the F/A conversion bv varying rhe pressure 
in the absorption cell. 

We looked tor signatures or chaotic be- 
havior with several numerical tools. We 
rtr^r examined the power spectrum (12). 
Periodic o>c illations in rhe data would pro- 
duce spikes in the power spectrum. The 
power spectrum from our dara i> v ery broad, 
indicating that the emissions are not purely 
periodic. The largest Lyapunov expo- 
nent — a measure of the rare at which tra- 
jectories initially close together in phase 
space diverge (12. 13) — is 0.2 3 > r 0.022 
for the filtered data in Fig. 1A, and the 
values ranged from 0.269 n 0.01 3 to 
0.314 — 0.01 1 in other measurements. Pos- 
itive Lyapunov exponents characterize cha- 
otic data, whereas periodic data have neg- 
ative values. There are two tests that mea- 
sure the dimension of the chaotic attractor 
in phase space: the capacity (or Hausdorff) 
dimension (12) and the correlation dimen- 
sion (14). Both measures yield a dimension 
or about 4 for the raw’ TDL data Laid about 



Fig. 2. A 3D time-delay phase pic: of the TDL 
frequency trajectory from the measurements dis* 
pia-.ea in Fig. 1A. The trajectory she-. vs intricate 
sc rai structure embedded within sp:m. structure. 
Ar animated version of this figure rev eo. s stretch- 
ing and folding motions. The points eoe colored to 
she .v the time evolution of the trajectcn from blue 
tc red. The axes indicate frequency m megahertz. 
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1.26S ± 0.070 for the filtered data). We 
also performed a single-value decomposi- 
tion (15) of our raw and filtered data, which 
resulted in only three significant eigenval- 
ues and eigenvectors, again indicating that 
the TDL trajectory has a low decree of 
complexity. 

A three-dimensional (3D) phase-space 
portrait (16) of the TPL filtered data using 
time-delay coordinates (Fig. 2) reveals a 
trajectory that is composed of strikingly 
complex embedded spirals. For stochastic 
data or random noise, such a plot shows no 
structure, and the trajectory randomly tills 
phase space. If the data are periodic, then 
the trajectory forms a closed orbit. An an- 
imated 3P phase portrait (17) shows the 
stretching and folding motions in the TPL 
frequency excursions. The TPL trajectory is 
confined, vet it never repeats — again a sig- 
nature of chaos. 

Because the TPL emission did appear 
chaotic, our next goal was to control this 
chaotic behavior. In 1990, Ott, Grebogi, 
and Yorke (OGY) (2) set the theoretical 
stage for chaos control based on the idea 
that the chaotic attractor of a system is 
made up of an infinite number of unstable 
orbits. The first step in the OGY method is 
to monitor one (or some) of the svstem 
variables and use the time-delay coordinate 
embedding technique (16) to construct the 
attractor for the system. From this attractor, 
they choose an orbit that yields the desired 
behavior and then locate the fixed points of 
that orbit. Finally, they apply a small time- 
dependent perturbation to an accessible sys- 
tem parameter to keep the fixed point in 
the path of the trajectory. The idea is to 
take advantage of the system s sensitive de- 
pendence on initial conditions and use mi- 
nor perturbations to control chaos. 
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Fig. 3. Diagram of the setup for our TDL chaos 
contro experiment. 


In 1 69 1 , Hunt developed a modification 
of the OGY method called occasional pro- 
portional feedback (OPF) to find, stabilize, 
and sw itch between different orbits of a di- 
ode resonator (7). One advantage of the 
OPF method is that detailed knowledge of 
the dynamics is not necessary, nor is it re- 
quired to know initially where the fixed 
points ot the system lie. In fact, the OI F 
method can he used to allow the system to 
find stable orbits by itself. The basic method 
of OrF is to sample an accessible dynamical 
system variable, and if the value of the vari- 
able falls within a prescribed range or “win- 
dow,” a system parameter is modulated with 
an amplitude proportional to the difference 
between the value of the sampled variable 
and the center of the window. An explana- 
tion of why OPF works and an extension of 
this method using the duration of the feed- 
back control as an additional control param- 
eter that is effective for higher complexity 
systems is found in (18). The circuitry need- 
ed to implement the OPF method is ex- 



Fig. 4. (A) A 2D time-delay phase plot of the 
chaotic (black) TDL and OPF -controlled (red) tra- 
jectories from Fig. 1 . (Bi Histogram of TDL trajec- 
tories in (A) with 0.2-Mhz bins that shows the 
extent o* the frequency variations. The red dash 
plot is a Gaussian fit to the OPF-controlled TDL 
data, Toe FWHM of the Gaussian, 1 .68 MHz. is an 
upper limit to the stabilized TDL frequency. The 
chaotic T DL excursion is - 1 8 MHz. 


t rente I v simple and easy to construct and has 
been successfully applied to a wide range of 
physical applications (9-9). 

For our problem, we monitor the fre- 
quency excursion, with the F/A method 
described above, and use the injection cur- 
rent of the TDL as the accessible system 
parameter (Fig. 3). The output frequency of 
a TDL is coarse-tuned by the operating 
temperature in a close-cycle refrigerator and 
fine-tuned by the injection current. The 
F/A output, from a detector preamplifier, 
enters the OPF electronics through a vari- 
able-gain input. An offset is added to the 
input signal to center it in a window. This 
signal is fed to a window detector and a 
sample-and-hold amplifier. The F/A signal 
is checked to determine if it lies within the 
adjustable width of the window. If so, the 
sampled signal is ted through a variable- 
gain output amplifier as a correction to the 
TDL control module, which adds it to the 
TDL injection current. The sampling rate 
and output correction signal are synchro- 
nized by an external clock, in our case 
running at 1 MHz. 

Control of the TDL emission was es- 
tablished quickly in a simple and robust 
fashion as the gain of the correction signal 
was gradually increased: The chaotic TDL 
emission (Fig. 1A) can be contrasted to 
the stability acquired when OPF control is 
fully active (Fig. IB). The controlled TDL 
emission is stable in both frequency and 
amplitude. 

Our measurements are all dominated by 
detector noise (the colored plots in Fig. 1). 
Because it is primarily high frequency (Fig. 
1C), we felt it appropriate to apply digital 
low-pass filtering to our data to improve our 
measurement accuracy. The amount of fil- 
tering was determined hv the requirement 
to preserve the characteristic of the TDL 
chaotic emission without any distortion. 

Several parameters determine how well 
the OPF method works. For example. TDL 
control is dependent on clock speed. For 
speeds slower than 600 kHz, OPF control 
fails; otherwise, the exact OPF clock speed 
is not critical: our experiments were all 
conducted at 1 MHz. This restraint i- con- 
sistent with (JO. H), and measurements 
here that show the time scale tor TDL 
variations is pis. The sensitivity ot the 
F/A conversion is another important tactor. 
High sensitivity in the F/A conversion, that 
is, small values of AY/Av (for voltage V and 
frequency v), or the steepness of the line 
where the laser is tuned causes difficulty in 
control. The ethylene line we used has a 12 
Y/GHz slope (to convert from voltage vari- 
ation to frequency, we used 81 MHz/Y). 
The w idth of the OPF window, on the other 
hand, is not a sensitive parameter. Once 
stable TDL emission was established, OPF 
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ainrmi wus robust and continued indefi- 
nitely without further adjustment. 

The slight drift in frequency in Fiy. IB is 
attributable to a small and slow 100-H: 
ramp in the injection current from the TDL 
control module. This drift shows that OPF 
control can he maintained while vary my 
the center frequency of the TLX. We 
achieved similar results with a fixed Fabry - 
Perot era Ion as the F/A be tun my the TDL 
to one side of a low-finesse etalon fringe. 
We maintained OPF control while scan- 
ning the frequency over a range of — SO 
MH:. This range is of practical importance 
as it demonstrates that OPF can he used 
with a TDL to scan across a molecular 
absorption feature and to perform spectro- 
scopic measurements at higher frequency 
resolution than previously possible. 

We can only yive an upper limit to the 
frequency stability attained with the OPF 
method because the measurements are 
dominated by detector noise that can only 
he improved by signal processing (Fiy. 1). 
Another factor is that our F/A process in- 
separably convolves amplitude and frequen- 
cy variations. To calculate an upper limit, 
we project the trajectories in two dimen- 
sions (Fiy. 4A). A histogram of the number 
of points in each 0.2-MH: bin of the con- 
trolled and uncontrolled emission is then 
made (Fiy. 4B). From a Gaussian fit to the 
controlled distribution, we obtain a value 
for its full width at half maximum 
(FWHM), here being 1.68 MH:, which we 
can use as a strict upper limit. 

An equal partitioning of the intensity 
and frequency variations in quadrature, a 
reasonable assumption, would reduce the fre- 
quence limit by V 2, to 1.19 MH:. The 
FWHM of the uncontrolled TDL excursion, 
measured by hand, is —IS MH:. Therefore, 
the OPF method improves the TDL frequen- 
cy stability by at least a factor of 11 to 15. 
We can determine more accurately the fre- 
quency stability by heterodyning the TLX 
output with a stabilized lO-jim CO z laser. 

The improved stability, in both the fre- 
quency and the amplitude of the OPF-con- 
trolled TLX, will enable new applications. 
We have mentioned higher resolution spec- 
troscopy that scans a TDL. An OPF-con- 
t rolled TDL can also be used as a local 
oscillator tor heterodyne radiometers. The 
far-infrared (FIR) spectral regime is ex- 
tremely important, and the recent develop- 
ment of wide-bandwidth hot-electron bo- 
lometer mixers makes the need tor a suit- 
able FIR local oscillator immediate. There 
are numerous candidates for future OPF 
control, such as an optically pumped sub- 
millimeter laser that emits from 30 p.m to l 
mm and whose ourput varies in a seemingly 
random fashion because of nonlinear feed- 
back between the pump and lasant gas. 
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ABSTRACT 

Close observations of the emission from a lead-salt tunable diode laser (TDL) 
reveal that the observed instrumental linewidth is actually a temporal aver- 
age of many narrow (less than 0.5 megahertz) emission "modes." Though the in- 
stantaneous laser output width is less than 0.5 MHz, the working width of the 
laser is limited to several MHz by frequency jitter. Since the jitter frequency is a 
dynamic system variable which occurs within a bounded range but varies in 
what appears to be a random fashion, chaotic behavior is suggested. With this 
in mind we designed a procedure to monitor and measure these frequency fluc- 
tuations on our own lead-salt TDL system. We then analyzed the data and dis- 
covered that the fluctuations are indeed indicative of a chaotic process. Utiliz- 
ing the chaos control technique of occasional proportional feedback (OPF) de- 
veloped in 1991 by E. R. Hunt, we constructed an electronic OPF controller. 
With this controller we succeeded in decreasing the frequency variations by a t 
least a factor of twenty over the same laser emission without the controller. 


1. Introduction 

In the 1990s, investigations of nonlinear dynamics have moved from purely 
theoretical and computational studies to experiment. Some of these experi- 
ments show that chaotic behavior is a useful feature in itself and more im- 
portantly may be exploited to perform practical functions 1 . Following theo- 
retical studies that reveal methods that stabilize unstable orbits 2 ' 3 a set of strik- 
ing experiments ensued that demonstrated that chaos can be controlled in a 
surprisingly simple manner 4 ' 5 and in as diverse a range of problems as cardiac 
arrhythmia 6 , thermal convection loops 5 , electronic circuitry 7 , an optical sys- 
tem 8 , and solid-state lasers 9,10 . We report here on new data and new results 
which validate and extend the work first reported in 1996 13 . 
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The practical demonstration of chaos control led us to consider 
whether these techniques are applicable to tunable diode lasers (TDLs). Lead- 
salt TDLs are valuable tools for the study of molecular parameters of atmos- 
pheric constituents important in planetary atmospheres 2 * * * * * * * * 11 . TDLs are operated 
at temperatures below 100 K and their frequency output is tuned by changing 
the diode temperature or injection current (which effectively causes a change 
in the temperature of the laser junction). 

Despite the extremely successful use of TDLs for laboratory spectros- 
copy, the full potential of a TDL has not been realized. Careful characteriza- 
tion of a TDL shows that the intrinsic linewidth of its emission is very nar- 
row (< 0.5 MHz or ~10' 5 cm' 1 ), a factor of 10 better than the instrumental 
resolution achieved in laboratory applications. The observed instrumental 
linewidth is actually a temporal average of many narrow frequency emissions 
or "modes" 11 ' 12 . These instantaneous modes, which are relatively constant in 
intensity, change on a time scale on the order of a microsecond. This behav- 
ior is a natural consequence of the intrinsically nonlinear operation of a laser 
junction where potentially many microscopic facets or lasing regions exist. 
Minute changes in the junction - such as in temperature, dimension, or index 
of refraction - can quench one lasing mode while allowing another to reach 
threshold. Though the instantaneous laser output width has been deter- 
mined to be less than 0.5 MHz, the working width of the laser is limited to 
several MHz 12 by frequency jitter resulting from this nonlinear operation. 
Since the jitter frequency is a dynamic system variable which occurs within a 
bounded range but varies in what appears to be a random fashion, chaotic be- 
havior is suggested. The detailed dynamics of such a complicated system is 
probably not amenable to analytical description. 


2. Frequency To Amplitude Converter - F/A 

Our first challenge is to capture the frequency fluctuation information from 

the TDL emission. Our solution is to construct a frequency-to-amplitude con- 

verter (F/A) with a fast response time. We use one of two approaches: either 

an absorption cell containing a gaseous molecular sample or a Fabry-Perot eta- 

lon is placed in the path of the TDL's beam. By adjusting the TDL to lase on 

the steep side of either the gas line or Fabry-Perot fringe, changes in the TDL 

frequency are converted to changes in detected signal amplitude. In Figure 1 

we see an illustration of this technique. The region on the side of an interfer- 

ometer fringe, or a molecular absorption line, looks approximately linear and 
has a fairly steep slope. For the sake of clarity in the figure, we show only the 

center portion of the side of the line, scaled for our convenience. We begin by 
tuning the laser frequency so that it is on the center of the side of the line. W e 
then note the amplitude of the signal there (marked as Center Line in the fig- 
ure). If the frequency of the laser decreases, as at point A, we immediately de- 
tect a decrease in the amplitude of the signal. Similarly, if the there is a de- 
crease in the laser frequency as at point B, we detect a decrease in signal ampli- 



tude. Having previously described the use of an absorption line F/ A 13 , in this 
paper we will examine the results of using an etalon fringe F/A. 



FIGURE 1 Frequency to Amplitude conversion. Shown is the side of an absorption line or etalon 
fringe. We tune the frequency so that it is at the center of the side at point C and hold the fre- 
quency there. If the frequency of the laser decreases, as at point A, we immediately see a de- 
crease in the amplitude of the signal. If the frequency increases, as at point B, we see an in- 
crease in amplitude. 


3. Examining The Frequency Fluctuations For Chaotic Character 

The F/A utilized in this investigation is a 3 " solid Germanium etalon, and 
we were operating the TDL at 1029 cm' 1 . In Figure 2 is plotted the frequency 
variation of the TDL emission. In (A) is plotted the raw data, and in (B) is 
plotted the same data filtered to remove noise. The data was taken with an 8- 
bit analog-to-digital converter, with a sampling rate of 500 kHz. 


Raw Data Filtered Data 



FIGURE 2 The frequency variation of a lead-salt TDL. In (A) is plotted the raw data, and in 
(B) is plotted the same data filtered to remove noise. 




We looked for the signature of chaotic behavior with several numerical tools. 
The largest Lyapunov exponent 14 was found to be 0.202 ± 0.055 for this data. 
The capacity dimension 15 was calculated to be 1.205 ± 0.066 for the raw data, 
and 1.368 ± 0.075 for the filtered data. In Figure 3 is plotted the correlation di- 
mension 16 as a function of embedding dimension for the Raw and Filtered 
data. We see that in the raw data, the presence of experimental background 
noise makes the correlation difficult to calculate, and also yields an artificially 
high value for the dimension. The filtered data gives a smoother plot, and 
shows a plateau onset for the correlation dimension of 2.858 ± 0.272. 



FIGURE 3 The Correlation dimension for the raw and filtered data. 


Another indicator of chaotic character is the power spectrum 15 . In Figure 4 we 
have plotted the log of the power spectrum as a function of frequency for the 
raw data. The power spectrum of the filtered data looks nearly identical except 
that it is attenuated at high frequencies. We see that the data has a very broad 
power spectrum indicating that the laser emission is not purely periodic. 



In Figure 5 we have used time-delay coordinate embedding 17 to create a 
phase-space portrait of the filtered data. The plot uses only 10,000 of the 65,000 



data points so that we may see the trajectory more clearly. The plot reveals a 
trajectory that is composed of strikingly complex embedded spirals. The trajec- 
tory is confined to a region of phase-space yet it never repeats. This is again 
indicative of a chaotic process. A phase plot of the entire data set may be seen 
in Figure 8 (A). 



FIGURE 5 Time-delay coordinate phase-space portrait for the filtered data. The axes indi- 
cate frequency in MHz. 


4. Control Of Chaos 

Having established to our satisfaction that the TDL emission indicates a cha- 

otic process, our next goal was the control of this chaotic behavior. In 1990 the 

first papers emerged which focused on the control of chaos in physical sys- 

tems. The first paper published by Ott, Grebogi, and Yorke 2 set the stage for the 

chaos control papers which followed. The method outlined in the paper, re- 

ferred to by most as the OGY method, was based on the idea that the chaotic 

attractor for the system is made up of an infinite number of unstable periodic 

orbits. The first step of the OGY method is to monitor one (or some) of the 

system variables, and use the delay-coordinate embedding technique 4 * * * * * * * * * * * * 17 to con- 

struct the chaotic attractor in the phase space of the system. Next, from this 

attractor find one of the orbits which yields the desired system behavior and 

locate the fixed points of that orbit. And finally a small time-dependent per- 

turbation is made in an accessible system parameter to keep the fixed point in 
the path of the system trajectory. This forces the system to remain in the de- 
sired orbit. The idea is to take advantage of the system’s sensitive dependence 
on initial conditions and use minor perturbations to the system to produce 
the desired behavior. 




Closely following the OGY paper Ditto, Rauseo, and Spano 4 used the 
OGY method to stabilize a physical chaotic two state system and switch be- 
tween the two states at will. These papers were followed quickly by a plethora 
of chaos control studies in widely diverse areas from controlling the flow of 
water in a thermal convection loop 5 to stabilizing cardiac arrhythmias 6 . 


5. Occasional Proportional Feedback 

In 1991 E. R. Hunt developed a modification of the OGY method which used 
occasional proportional feedback (OPF) to find, stabilize, and switch between 
nineteen different orbits from the chaotic attractor of a diode resonator 7 * . One 
advantage of the OPF method is that we do not need a detailed knowledge of 
the dynamics of the system nor is it necessary to know initially where the 
fixed points of the system lie. In fact, the OPF method can be used to allow the 
system to find stable orbits by itself. The basic method of OPF is to sample an 
accessible dynamic system variable, and if the value of the variable falls 
within a prescribed range or "window", a system parameter is modulated 
with an amplitude proportional to the difference between the value of the 
sampled variable and the center of the window. 

Another difference between the OPF method and the OGY method is 
that the perturbations made to the accessible system parameter, necessary to 
stabilize the orbit, are not restricted to being small in the OPF method. Hunt 
found that with large perturbations the attractor could be altered and periodic 
orbits could be found where none had existed previously. OPF has more re- 
cently been used to control a chaotic solid state laser system 9 and also to syn- 
chronize a pair of diode resonators 18 . An explanation of why OPF works and 
an extension of this method using the duration of the feedback control as an 
additional control parameter that is effective for higher complexity systems is 
found in Carr and Schwartz 19 . 


6. Our OPF Controller System 

For our OPF controller, we monitor the frequency excursion of the laser out- 

put (with the F/A method described above) as our dynamic system variable, 

and use the injection current to the TDL as our accessible system parameter. 
In Figure 6 is a diagram for our OPF controller and TDL system. The output 
frequency of the TDL is coarse-tuned by the operating temperature in a close- 
cycle helium refrigerator and fine-tuned by the injection current. The F/A 
output, from a detector pre-amplifier, enters the OPF electronics through a 
variable-gain input. An offset is added to the input signal to center it in a 
window. This signal is fed to a window detector and a sample-and-hold am- 
plifier. The F/A signal is checked to determine if it lies within the adjustable 
width of the window. If so, the sampled signal is fed through a variable-gain 
output amplifier as a correction to the laser current supply where it is added 
to the TDL injection current. The sampling rate of the OPF controller sample- 



and-hold amplifier and the timing of the output correction signal are syn- 
chronized by an external clock, in our case running at 1 MHz. 



FIGURE 6 OPF controller and system diagram. 


Several parameters determine how well the OPF method works. For 
example, TDL control is dependent on clock speed. For clock speeds slower 
than 600 kHz, OPF control diminishes. Otherwise the exact OPF clock speed is 
not critical and our experiments are all conducted at 1 MHz. The sensitivity of 
the F/ A conversion is another important factor. High sensitivity in the F/A 

conversion, that is small values of AV/Av (for voltage V and frequency n), or 
the steepness of the line where the laser is tuned, causes difficulty in control 13 . 
The etalon fringe we used here has a 6.7 V/GHz slope (to convert from volt- 
age variation to frequency, we used 150 MHz/V). This slope is much lower 
than the slope of the molecular absorption line used in previous experi- 
ments 13 , and the control is smoother. The width of the OPF window, on the 
other hand, is not a sensitive parameter. Once stable TDL emission was estab- 
lished, OPF control was robust and continued indefinitely without requiring 
any further adjustments. 




7. 


Experimental Results 


Figure 7 shows the results from one of our experiments. Control of the TDL 
emission was established quickly in a simple and robust fashion as the gain of 
the correction signal is gradually increased. The chaotic TDL emission (Fig. 7 
A) can be contrasted to the stability acquired when OPF control is fully active 
(Fig. 7 B). (The uncontrolled data (A) is the same data as in Figure 2). The con- 
trolled TDL emission is stable in both frequency and amplitude. 



data was filtered to remove the detector noise. The interval of each measurement was 130 msec, 
and the data in the three panels were not taken contiguously. 


Figure 7 C shows the detector noise, here the TDL laser beam is blocked. Our 
measurements are all dominated by detector noise. Because the detector noise 
is primarily high frequency, we felt it appropriate to apply digital low-pass fil- 
tering to our data to improve our measurement accuracy. Therefore the data 
shown in Figure 7 shows the results of the filtering. (The raw uncontrolled 
emission can be seen in Figure 2.) The amount of filtering was determined by 
the requirement to preserve the characteristic of the TDL chaotic emission 
without any distortion. 

We shall now compare the frequency fluctuations in the controlled 
versus the uncontrolled data. The total deviation for the uncontrolled laser 
emission is 109 MHz, with a standard deviation of 17.3 MHz for the unfiltered 
data, and 104 MHz, with a standard deviation of 17.2 MHz for the filtered data. 
We see that the filtering of the high frequency noise has very little effect on 
the uncontrolled laser emission (as can be seen visually in Figure 2). 

For the controlled laser we find that the total deviation is 11.7 MHz, 
with a standard deviation of 1.45 MHz for the unfiltered data, and 5.57 MHz, 
with a standard deviation of 0.63 MHz for the filtered data. While the post-fil- 
tering of the data had very little effect on the uncontrolled laser emission, we 
see that the standard deviation of the controlled emission is reduced by over 
half, indicating that the controlled laser signal is background noise limited. 



The filtered background signal noise has a total deviation of 11.7 MHz, 
with a standard deviation of 1.35 MHz for the unfiltered data, and 6.08 MHz, 
with a standard deviation of 0.83 MHz for the filtered data. The comparison of 
the above statistics for the controlled emission and the background detector 
noise again indicates that the controlled laser emission is background noise 
limited, and that our filtering is therefore justified. For the filtered data, the 
ratio of the standard deviations for the uncontrolled and controlled emission 
indicate a reduction in frequency fluctuation by a factor of 27. 


8. Upper Limit To Frequency Stabilization 

We can only give an upper limit to the frequency stability attained with the 
OPF method because the measurements are dominated by detector noise that 
can only be improved by signal processing. Another factor is that our F/A 
process inseparably convolves amplitude and frequency variations. To calcu- 
late an upper limit, in Figure 8 (A) we project the trajectories in two dimen- 
sions of the uncontrolled (black) and stabilized (small white region in center) 
laser emission. In Figure 8 (B) we have plotted the histograms for both the 
uncontrolled laser emission data (1), and the controlled emission data (2). The 
bin width for the histogram is 1 MHz. A Gaussian profile was then fitted to 
each histogram, and the profile was over-plotted onto the histogram. The 
Gaussian fitted to the uncontrolled laser data has a FWHM of 29.5 MHz, 
while the Gaussian fitted to the controlled laser data has a FWHM of 1.45 
MHz. We may use this FWHM of 1.45 MHz as a strict upper limit. An equal 
partitioning of the intensity and frequency variations in quadrature, a reason- 
able assumption, would reduce the frequency limit by ^2 or to 1.03 MHz. 
Since the FWHM of the uncontrolled TDL excursion was found above to be 
29.5 MHz, we find that the OPF method improves the TDL frequency stability 
by at least a factor of 20 to 29. 


9. Scanning The Center Frequency Of The Controlled TDL 

In Figure 9 we have plotted the successful scan of the center frequency of the 
TDL with the OPF controller activated. The range of the scan is 120 MHz (4 
mem' 1 ). This is of practical importance as it demonstrates that OPF can be used 
with a TDL to scan across a molecular absorption feature and to perform spec- 
troscopic measurements at higher frequency resolution than previously pos- 
sible. 

The improved stability, both in frequency and amplitude of the OPF- 
controlled TDL, will enable new applications. In addition to the high resolu- 
tion spectroscopy mentioned above, an OPF-controlled TDL can also be used 
as a local oscillator (LO) for heterodyne radiometers. The far-infrared (FIR) 
spectral regime is extremely important and the recent development of wide- 
bandwidth hot-electron bolometer mixers makes the need for suitable FIR LO 





FIGURE 8 In (A) is a 2D time-delay phase plot of the chaotic uncontrolled (black) TDL and 
the OPF controlled (small white region at center) trajectories from the data in Figure 7. In (B) 
are histograms of the TDL trajectories in (A) (uncontrolled (1), controlled (2) ) with bin widths 
of 1 MHz. Each histogram is then over-plotted with a Gaussian profile fit to the histogram 
data. For the uncontrolled emission the FWHM of the fitted Gaussian is 29.5 MHz, while for 
the controlled emission the FWHM is 1.45 MHz. This width of 1.45 MHz is an upper limit to 
the stabilized TDL frequency. The Ratio of the uncontrolled and controlled FWHM indicate an 
effective narrowing of the laser frequency output by a factor of 20. 





immediate. There are numerous candidates for future OPF control, such as an 
optically pumped submillimeter laser that emits from 30 p.m to 1 mm and 
whose output varies in a seemingly random fashion because of nonlinear 
feedback between the pump and lasant gas. 



controller activated. The ramp generator used was set for a 10 msec (100 Hz) scan, where the 
scan stays at the minimum value for 5 msec then increases for 5 msec before returning to the 
minimum value. At the right edge we can see where the OPF controller is retarding what would 
otherwise be a nearly vertical fall of the frequency back to the minimum value. 
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